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Abstract. We study the dynamics of a classical, two-component plasma in two dimensions, in the 
vicinity of the Kosterlitz-Thouless (KT) transition where the system passes from a dielectric low- 
temperature phase (consisting of bound pairs) to a conducting phase. We use two "complementary" 
analytical approaches and compare to simulations. The conventional, "intuitive" approach is built on 
the KT picture of independently relaxing, bound pairs. A more formal approach, working with Mori 
projected dynamic correlation functions, avoids to assume the pair picture from the start. We discuss 
successes and failures of both approaches, and suggest a way to combine the advantages of both. 



1. INTRODUCTION 



*— > 

c^ , The two-component plasma in two dimensions (2D) is the generic model for topological excitations in 
2D systems with a U(l) order parameter symmetry, usually called vortices |1|. Prominent examples 
of such systems are the "2D superfiuids" , superconducting or ^He films and Josephson junction 
arrays (JJA's), others are the XY model of 2D planar magnets, and — with some modifications — 
2D melting 0. The 2D plasma undergoes the well-known Kosterlitz-Thouless (KT) transition at a 
finite temperature Tkt- 

The dynamic behavior close to this transition is of both principal and practical interest. It is 
conveniently described by a frequency dependent, complex dielectric function e{uj), related to the 
dynamic correlation function ^p{k,uj) of the charge density p{r,t) of the plasma: 

q^nSp{k) + ico'$p(fc, uj) 



oo 

%ik,u) = |dte'-*|dVe-''=-^(p(r,t)p(O,0)) . (2) 



,^ . Here n is the mean number density of the particles with charge ±g, and Sp{k) = $p(fc,t=0)/g n 
^ '. is the static charge structure factor, determining the zero frequency limit of e{uj). As usual, e^u) is 
related to the dynamic conductivity criuj) by e{uj) = 1 — 27icr{uj)/iu. 

Recent numerical simulations by Jonsson and Minnhagen (JM) [Q of the 2D XY model with 
Ginzburg-Landau dynamics have shed new light onto the critical dynamics of the associated vortex 
plasma. Using an appropriate definition for the vortex density p{r), JM have computed <l>p and 
extracted the complex dielectric function e{uj) using ([I|). It shows the following, interesting features: 

— l/e{uj) approximately has a scaling form with a characteristic frequency cUch going to zero as 
T -^ Tkt ( "critical slowing down" ) from both sides. 

— For low frequencies, Re[l/e(ct;) — l/e(0)] oc u over a relatively large ut interval, contrary to a 
simple "Drude" behavior (oc u"^). This corresponds to a non-analyticity (j{uj) oc ln\uj\ at low u. 

JM obtain very good fits of their data by a phenomenological expression proposed earlier |I| by 
Minnhagen. Experiments on JJA's 0] seem indeed to confirm the non-analytic frequency dependence 
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Figure 1: Left: Real and imaginary parts of e{£,)/e{Lo) plotted against a scaled frequency for a series of 
different temperatures above the KT transition (full curves, to = —0.18, —0.15, —0.09, 0,0.09, and ^/rQ = 
6 X 10^^,9134,63,8.7,3.7 from above to below). The scaling function is almost indistinguishable from the 
simple phenomenology of JM (dashed). For comparison, a simple Drude form is also shown (dotted). 
Right: The characteristic frequency uo^h defined in the text as a function of temperature T, above as well as 
below Tkt (= vertical dashed line). 



of e{uj), and various theoretical attempts have been made 0] in order to explain this finding. The 
critical slowing down has, to our knowledge, not been theoretically addressed so far. 



2. AHNS PAIR RELAXATION PICTURE 

The "standard" dynamic theory of the KT transition has been developed in 1978 by Ambegaokar 
et al. (AHNS) 0, who built up the dielectric function e{uj) additively from contributions of pairs of 
any size and of free particles. We use a slight modification of this idea here, defining by 
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a dielectric function of bound pairs, and setting e{uj) = e^iui) — 2nq'^ni/iujT. The scale dependent 
dielectric constant e(r) is determined by the KT flow equations [|l], 0, and the factor 1/[1 — iu;r(r)] 
in the integral describes the relaxation of individual pairs of size r. ro is a minimum length of 
the order of the vortex core size, ^ is the KT screening length, and Uf ~ ^^^ the density of free 
particles. To determine e(r), we use Minnhagen's fl]] version of the KT flow equations for the scale 
dependent "reduced temperature" t(r) = 1 — q^ /Al{r)kjiT and fugacity y{r). Their solutions are 
known analytically, in the regions T < Tkt and T > Tkt separately [H. The integral (^ is done 
numerically. 

We flnd that the frequency at which Im[l/e(ci;)] has its maximum is more or less temperature 
independent below Tkt, whereas above it varies rapidly like ^(T)^^. Both above and sufficiently far 
below Tkt, the low-o; behavior of Re[l/e] is closer to linear than a simple Drude form. However, 
above Tkt there flnally is a crossover to "Drude" uj"^ behavior as cj — >■ 0, and below Tkt both real and 
imaginary part of l/e(u;) — l/e(0) behave as cj" at low uo, with an exponent a = 2e(n)k t ~ ^ which 
goes to zero as T — > Tkt- If normalized by a factor of e(0) our data above Tkt display an almost 
perfect scaling behavior with a characteristic frequency 0Jch{C) ~ 10Q'^/^^^(Or, over an extremely 
wide range in ^ (Fig. ||L). We find no comparable scaling below Tkt, but if cJch is defined here by the 
condition Re[l/e(to'ch) — l/e(0)] = — | Im[l/e(u;ch)] (as done by JM), it decreases rapidly as T — ^ Tkt 
from both sides. In Fig. ^R we have plotted the corresponding temperature dependences of cjch both 
above and below Tkt, which is in good qualitative agreement with Fig. 6 of JM [Q. 
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Figure 2: Left: Real and imaginary parts of the reduced charge friction function 7(ci)), calculated from Eq. 
(0) with A = 2. Right: The resulting l/e(a;) (fuh curves) compared to the phenomenology of JM (dashed) 
and to a Drude Form (dotted), with maxima of — Im at the same frequency {uj ~ 0.22) and amplitudes 
adjusted. 



3. CALCULATIONS USING MORI'S TECHNIQUE 

Mori's technique [0] for evaluating dynamic correlation functions yields for the charge density corre- 
lator at small wave numbers k: 



^p{k,u) 



q^nSp{k) 



—lUJ + 



kjjTk'^ 



(4) 



5p(fc)7p(a;) 

The "memory function" (or "random force" correlation function) 7^(0;) obeys 
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(5) 



where Sn^it) denotes the deviation of the local number density from its mean value, and V is the 
system volume. 7p(co') generalizes the bare friction constant F, adding a contribution due to the 
particle interaction. 

In the spirit of current practice in the theory of liquid dynamics [0] we factorize, as a first step, 
the combined correlator into a product of a number and a charge correlation function: 



{Sn_k{t)pk{t)6ni{0)p_i{0)) ^ V^6M'^n{-k,t)<^p{k,t) 



(6) 



Using ^n{k,uj) ^ n/{—iuj + MZ^) and Sp{k) ^ k'^/{k'^ + fe ^) (the length 6 is a measure of the 
interparticle distance) we obtain the implicit equation 
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(7) 



determining the function 7(ci)) = 7p(u;)/F which depends on cD = cj/cUch,!- The characteristic fre- 



quency ujch,i = k^T/b T reflects the spatial density of particles, and A = Z^ ™.2 is a dimensionless 
parameter of order unity. 

Fig. I^L shows real and imaginary parts of the solution of (0) above Tkt, where cr(u;) oc ■jp{uj)~^. 
Re7 shows the expected (see Sec. 1) logarithmic uj dependence in some region, but has a finite limit 
for tu — > 0. As seen in Fig. ^R, our mode decoupUng approach yields a dielectric function somewhat 
closer to the phenomenology of JM than a pure Drude form. However, the critical slowing down is 



missing, since the factorization approximation (P) underestimates effects of charge binding into pairs 
on length scales < ^. 

A simple — and admittedly ad hoc — way to build pairing into the dynamic friction function 
consists in adding to the r.h.s. of (|^) a slowly decaying term oc e~*/'^°='', corresponding to a contribution 
Tpairsli^) = 7o/[Te^c ~ ^^1 ^^ 7p(^) ^^^ describing the random force correlations of pairs that unbind 
after a typical time Tesc- We interprete Xesc as the "thermal escape time" out of the logarithmic pair 
potential, screened at a distance ^: 

resc(0«^f-y ' - = 7^- (8) 

q^ \roJ KT 

The dynamic exponent z ~ 4.7 in the vicinity of the KT transition. Studying the correlations in 
an ensemble of independent pairs with average squared size (av^) pa 46^, the amplitude of 7pairs 
can be estimated as 70 ~ q^ / k^T {at'^) . The new term 7pairs dominates 7^ for large ^ and low u;, 
and so the static limit of the solution diverges as 7p(c<j=0) oc resc(0 0^ ^^- As a consequence, both 
the new characteristic frequency c<Jch,2 = ^e^c ^^^ the static conductivity o"(0) = q^n/'jpiO) of the 
plasma vanish smoothly as ^(T)~^. The prediction 2; ?^ 4.7 is in striking contrast to the scaling 
with 2; = 2 in Fig. |1|L obtained on the basis of AHNS theory, and to the usual expectation that 
(t(0) ~ q^Uf/T ~ ^~^. We note, however, that recent dynamical scaling analyses [§ for different 
realizations of 2D superfluids actually do suggest anomalously large dynamical exponents in the 
range z k, 4.5 ... 5.9. 

4. CONCLUSIONS 

We have studied here the critical dynamics of a 2D, two-component plasma in the vicinity of the 
KT phase transition, using a modification of AHNS theory and a Mori type approach. The AHNS 
approach is quite successful as a convenient phenomenological interpretation of simulation data, 
but conceptually dissatisfactory in that pairs are put in "by hand". The Mori approach is more 
fundamental, starting from microscopic equations of motion, but the correlations it predicts (within 
a mode decoupling approximation) are not strong enough to imply critical slowing down. We have 
proposed a way to cure this deficiency by introducing additional, long-lived random force correlations, 
describing temporary pairing. A virtue of this extension is that it contains two frequency scales, one 
finite and the other vanishing at Tkt, in agreement with the numerical findings. 
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